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INTRODUCTION 


A general formulation for the analysis of steady and un- 
steady, subsonic and supersonic potential aerodynamics for 
arbitrary complex geometries Is presented here. The proposed 
paper Includes the theoretical formulation, the numerical pro- 
cedure and numerical results. In particular, generalized 
forces for fully unsteady (complex frequency) aerodynamics 
for an AGARD coplanar wing-tall Interfering configuration In 
both subsonic and supersonic flows are Included In the paper. 

The theoretical formulation Is based upon an Integral 
equation presented In Refs. 1 and 2, which Includes completely 
arbitrary motion. Steady and oscillatory aerodynamic flows 
are considered In Refs. 3 and 4 (enclosed here). A review of 
the problem Is given In Ref. 4 and therefore Is not Included 
here. 


Here a much more general formulation Is considered. 
First, small-amplltude, fully transient response In the time 
domain Is considered. This yields the aerodynamic transfer 
function (Laplace transform of the fully unsteady operator) 
for frequency domain analysis (Ref. 5 enclosed here). This 
is particularly convenient for the linear systems analysis 
of the whole aircraft. The formulation briefly outlined in 
Ref. 5 has now been completed and implemented in the computer 
program SOUSSA (Ref. 6, for subsonic and supersonic). 
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In addition, no diaphragm Is required In supersonic flow, 
so that the subsonic and supersonic formulations have been 
unified. Also the formulation has been extended to allow for 
the analysis of the problem of xne wing-wake Intersecting the 
tall. Additional features are described In details In the 
section entitled "Assessment of 1/ethod". 

The new formulation, program and results will be fully 
described In the proposed paper. 


METHOD OF SOLUTION 

The method presented here is based upon a formulation 
developed by Morino.^ For simplicity, only the incompres- 
sible steady state is briefly described here. The formula- 
tion, by making use of the Green function method applied to 
the equation of the velocity potential, yields an integral 
equation relating the unknown potential on the surface of 
the body to its known normalwash. By making use of the 
finite-element method, and by the assumption that the poten- 
tial is constant within each quadrilateral element, the in- 
tegral equation is approximated by a linear system of N 
equations relating N (unknown) values of the potential to N 
(known) values of normalwash at the centroids of N elements. 

For the sake of generality and flexibility, in particu- 
lar, for structural analysis, the normalwash is expressed in 
terms of the generalized coordinates and generalized velocities 


From the potentials at centroids of elements, by an 
averaging scheme (by which the potential at a corner Is 
approximated by the average value of potentials at the 
centroids of the elements In Its Immediate surroundings), 
the potentials at the nodal points are obtained and hence the 
potential at any point on the surface can be expressed by a 
finite-element Interpolating formulation with bl-llnear local 
shape functions. Finally, the pressure coefficients and 
generalized forces can be evaluated by a simple finite- 
element procedure. 


NUMERICAL RESULTS 


Typical numerical results obtained with SOUSSA are pre- 
sented in this section. Due to space limitations, the result 
are only very briefly outlined. 

Figures 1 and 2 are the l^ft and moment coefficients 
of a rectangular wing oscillating in pitch with 
Mach number ranging from 0 to 2.5. Results for 
the supersonic flow were obtained without the use 
of diaphragms and have never been presented before. 

The comparison against Ref. 11 is in general, in 
excellent agreements. Figures 3, 4 and 5 present 
the pressure distributions of a rectangular wing 
in steady subsonic and supersonic flow, and again 
the results are in very good agreemert with the 
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on«8 of Ref. 9. Figures 6, 7 and 8 are results for 
a wing-body configuration in both steady and fully 
unsteady flow, for both subsonic and supersonic 
speeds. 

Figures 6 and 7 are presented in order to demon- 
strate the capability of the present method of 
analyzing fully unsteady flows. Figures 9 and 10 
Include the results for simple wings with control 
surface in steady and oscillatory flows. Figure 11 
presents flutter applications ( in excellent agree- 
ment with the results of Ref. 17). Tables 1 through 
3 are the generalized forces for an AGARD wing-tail 
configuration in quasi-steady and oscillatory flow 5n 
comparison with existing methods. 

Further results, such as the fully unsteady aero- 
dynamic analysis of the AGARD wing-tall configuration 
and other wing-body-tail configurations (with control 
surfaces) will be included in the proposed paper. 


ASSESSMENT OF METHOD 


An assessment of the method is briefly outlined here. In 
particular, unique features of the methodology (with existing 
methods) are stressed and progress with Ref. 4 is emphasized. 
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It provides a uzxlfied approach for steady, oscilla- 
tory and fully unsteady subsonic and supersonic 
aerodynamic floes* 

To the authors' knowledge, this Is the only exist- 
ing program which can handle fully unsteady (complex 
frequency) aerodynamics for complex configurations 
(e.g., wlng-hody-tall combination) In both subsonic 
and supersonic regimes. No other program can handle 
even oscillatory supersonic flow problems for 
complex configurations. 

It can be applied to arbitrarily-complex configur- 
ations. Wing -body -tail configurations with control 
surface have been analyzed. (No existing results 
are available for comparison. However, results 
for a simple wing with control surface shows that 
the present method Is In good agreement with 
existing ones.) 

It Is a unique feature of the present method In 
that it Is capable of analyzing influence of 
actual wake geometry on aircraft surfaces. In 
particular, it can analyze wing-tail configura- 
tions with the wake generated by the wing inter- 
cepting the tall and merging with the tall gen- 
erated wake when it leaves the tall trailing edge. 
This feature has been added very recently to the 
program and Improves considerably the range of 
applicability of the code. Results are in good 
agreement with existing ones. 


6 . 


$. The program ia computationally extremely general, 
flexible, efficient and above all, accurate. The 
elimination of diaphragms in supersonic flow improved 
considerably the simplicity and efficiency of the 
code. (Ref. 4 requires the use of diaphragm and 
hence is limited to simple geometries). 

6. In contrast to existing methods, which in many in- 
stances require extensive user's background in 
aerodynamics and familiarity with the specific 
method, the present code requires limited human 
Itervention and is very easy to use. 

7. The evaluation of pressure is performed using the 
finite-element method. (Ref. 4 used finite- 
difference and was limited to thin wings) 

8. The generalized forces are evaluated for arbitrary 
geometry and arbitrary three dimensional mode shapes. 

9. Another unique feature of the present method on un- 

steady potential flow problems in that the flutter 
analysis often requires the analysis on a specific 
geometry for a wide range of frequencies. In the 
present method, the frequency-dependent coefficients 
of the aerodynamic transfer matrix, is expressed as 
a combination of complex frequency-independent 
coefficients* with simple frequency-dependent 
coefficients: the advantage is that every addi- 

tional frequency analysis other than the first 

one requires only a minimal amount of CPU time 




coefficients nf Ref. enclosed here 
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10. In lt«r«tlve proosdurts (for Instance for optimal 
design) It Is generally required to predict 
generalized aerodynamic loads due to a variety of 
vibration modes. In tbe present method, the 
aerodynamic coefficient matrix Is written as the 
product of three matrices. The first and the third 
(for the normalwash and for the evaluation of the 
generalized forces) are mode dependent but very 
simple, while the second one (relating pressure 
distribution to normalwash distribution) Is mode 
Independent. By the same reasoning as above, 

the CPU time required for additional modal analysis 
Is reduced to a relatively negligible level. 

11. Evaluation of the normalwash for complex configura- 
tions from prescribed three dimensional mode shapes 
(Ref. 4 was limited to thin wings with vertical 
displacements.) Is available. Oownwssh due to 
turbulence Is also included. 

12. Applications to flutter have been considered. The 
results (see next section) are in good agreement 
with existing ones. 

In conclusion, the proposed paper will emphasize the 
generality, flexibility, efficiency of the present method. 
Last, but not least, tbe present method provides a unified 
approach to cover the whole linearized potential flow spectrum 
with very limited human Intervention required in using the 
computer code SOUSSA. 
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— — — two dimensional airroil 
0 present method 


Figure / Flutter speed as a Function of ® for a rectangular wing with 
AR ~ 16, M — 0, T = 0, 1%, 4 = 5, Xof = 0.2, — 0.5, and 

NX = 8 , NY = 10, Results are compared with exact solution 
given by tvro dimensional airfoil theory (Ref jlo.) O^ca = -0.2C] 
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0 present method 


Figure f/b Flutter speed os o function of ° rectangular wing with 

AR = 16, M = 0, T = 0.1%, u = 10, Xa ^ 0.2, r» = 0.5 and 
NX = 8 , NY = 10. Results ore compa ed with exoct solution 
given by Kvo dimensional airfoil theory (Ref.Jo)(XcA = -0.2C), 
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